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Abstract. The efficiency of the future devices for quantum information processing 
is limited mostly by the finite decoherence rates of the qubits. Recently a substantial 
progress was achieved in enhancing the time, which a solid-state qubit demonstrates 
a coherent dynamics. This progress is based mostly on a successful isolation of the 
qubits from external decoherence sources. Under these conditions the material-inherent 
sources of noise start to play a crucial role. In most cases the noise that quantum device 
demonstrate has 1// spectrum. This suggests that the environment that destroys the 
phase coherence of the qubit can be thought of as a system of two-state fluctuators, 
which experience random hops between their states. In this short review we discuss 
the current state of the theory of the decoherence due to the qubit interaction with 
the fluctuators. We describe the effect of such an environment on different protocols 
of the qubit manipulations - free induction and echo signal. It turns out that in many 
important cases the noise produced by the fluctuators is non-Gaussian. Consequently 
the results of the interaction of the qubit with the fluctuators are not determined by 
the pair correlation function only. 

We describe the effect of the fluctuators using so-called spin-fluctuator model. Being 
quite realistic this model allows one to evaluate the qubit dynamics in the presence of 
one fiuctuator exactly. This solution is found, and its features, including non-Gaussian 
effects are analyzed in details. We extend this consideration for the systems of large 
number of fluctuators, which interact with the qubit and lead to the 1// noise. We 
discuss existing experiments on the Josephson qubit manipulation and try to identify 
non-Gaussian behavior. 
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1. Introduction 

Coherence in quantum solid state devices inevitably suffers from fluctuations due to their 
environment. In particular, rearrangement of electrons between traps in the insulating 
regions of the device, as well as stray flux tubes causes pronounced fluctuations in many 
quantum devices. At low frequencies part of these fluctuations typically has a 1/f 
spectrum and is referred to as 1// noise. Such noise is generic for all disordered materials 
(for a review see, e.g., [H E]), it is also common in single-electron and other tunneling 
devices, see, e.g., [3]. Recent experiments [H [5] on Josephson qubits indicated that 
charged impurities may also be responsible for noise. Low frequency noise is specifically 
harmful since it is difficult to filter it out by finite band filters. 

One of the most common sources of low-frequency noise is the rearrangement of 
dynamic two-state defects, fluctuators, see, e. g., the book [6] and references therein. 
Random switching of a fluctuator between its two metastable states (1 and 2) produces 
random telegraph noise. The process is characterized by the switching rates 712 and 
721 for the transitions 1 — * 2 and 2^1. Only the fluctuators with energy splitting E 
less than temperature, T, contribute to the dephasing since the fluctuators with large 
level splitting are frozen in their ground states (we measure temperature in the energy 
units). As long as E <T the rates 712 and 721 are close in magnitude, and without loss 
of generality one can assume that 712 ~ 721 = 7- i- e., the fluctuations can be described 
as a random telegraph process (RTF), for reviews see [SI [7]. A set of random telegraph 
fluctuators with exponentially broad distribution of relaxation rates, 7, produces noise 
with 1/f power spectrum at 7min <^ uj = 2-nf <ti 70. Here 7min is the switching rate of 
the "slowest" fluctuator while 70 is the maximal switching rate for the energy difference 
between the fluctuators's metastable states equal to temperature. Random telegraph 
noise has been observed in numerous nanodevices based both on semiconductors, normal 
metals, and superconductors [8]. 

Specific features observed in recent experiments [9] in Josephson phase qubits were 
interpreted in terms of resonant interaction of the qubit with two-level impurities. 
These experiments, as well as results [10] in a superconducting quantum circuit 
(quantronium), stimulated a very important attempt [11] to establish a relation between 
the contributions of two-level fluctuators with low and high frequencies. That has 
become possible because both of the frequency domains contribute to the dephasing time 
and the energy relaxation time, respectively. However, this relation was significantly 
based on the assumption that the statistics of the fluctuations of the qubit parameters 
are Gaussian. In our view, this assumption is not obvious, and the main aim of this 
short review is to discuss he applicability range of the Gaussian approximation, as well 
as the deviations from the Gaussian approximation in connection with the problem of 
qubit dephasing. For this purpose, we will use a simple classical model within which one 
can evaluate exactly the qubit response to typical manipulation protocols. This model 
if often referred to as the spin-fluctuator (SF) model in similarity with the widely used 
and sometimes overused spin-boson model. According to the SF model, the quantum 
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system - qubit - interacts with a set of two-level entities. The latter stochastically 
fluctuate between their states due to interaction with a thermal bath, which may be 
not directly coupled with the qubit. Since we are interested in the low-frequency noise 
generated by these switches the latter can be considered as classical. Consequently, the 
system qubit -|-fiuctuators can be described by relatively simple stochastic differential 
equations, which in many cases can be exactly solved. In particular, many results can 
be just borrowed from much earlier papers on magnetic resonance [T2l [13], on spectral 
diffusion in glasses [T4j, as well as works on single- molecule spectroscopy [T5] . 

The SF model has previously been used for the description of effects of noise in 
various systems [T6],[I71[T8],ll9ll2ni[2l] and was recently applied to analysis of decoherence 
in charge qubits [22[23l[25l[2l[26l[27l[28l[29l[30]. Various quantum and non-Markovian 
aspects of the model were addressed in Ref. [31]. These studies demonstrated, in 
particular, that the SF model is suitable for the study of non-Gaussian effects and 
that the latter may be essential in certain situations. 

In this short review we will address non-Gaussian effects in decoherence of 
Josephson qubits using the echo signal as an example. Resonant interaction between the 
qubit and fluctuators discussed in P [321 [33] will be left out. The paper is organized as 
follows. In Sec. [2] we briefly describe decoherence in the Gaussian approximation. The 
spin-fiuctuator model is described in Sec. [31 Section [His aimed at discussion of relevance 
of the SF model to existing experiments, while current microscopic understanding of the 
two-level fluctuators is briefly reviewed in Sec. [51 



2. Gaussian decoherence 

A qubit is described by the generic Hamiltonian of a pseudospin 1/2 in a "magnetic 
field" B, which can be time-dependent: 

n, = lB-a (1) 

where cTj are the Pauli matrices. It is well known that any state vector, |\E'), of the qubit 
determines the Bloch vector M through the density matrix 

p = |M/)(M/| = i(l + M-^). (2) 

The Schrodinger equation turns out to be equivalent to the precession equation for the 
Bloch vector: 

M = B X M. (3) 

The problem of decoherence arises when the "magnetic field" is a sum of a controlled 
part Bq and a fluctuating part h{t) which represents the noise, i.e., is a stochastic 
process determined by its statistical properties, B = Bq + b(t). The controlled part Bq 
is not purely static - to manipulate the qubit one has to apply certain high-frequency 
pulses of Bo in addition to the static fields applied between manipulation steps. In this 
paper we will always assume that the manipulation pulses are short enough and neglect 
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Figure 1. Bloch vector in the rotating frame of reference. 



the decoherence during the pulses. Therefore we need to consider only the effect of noise 
in the presence of static Bq. 

We will consider only the case where b || Bq and let the 2;-axis lie along the common 
direction of Bq and b, see Fig. [1] This situation is called pure dephasing because z- 
component of the Bloch vector, M^, is conserved during the process. As long as the time 
evolution of M is governed by Eq. ([3]), the length |M| is also conserved, while the length 
I (M) I of the vector M averaged over the stochastic process b decays. Description of this 
decay is the main objective of the decoherence theory. In the case of pure dephasing, 
this will be the decay of the components and My. It is convenient to introduce a 
complex combination = {M^ + iMy)/ y/MlTMl- Equation ^ can be written in 
terms of m_|_ as 

?7i+ = iBm+ with solution m+(t) = e-^o '^^*'^'^*'m+(0). 

This solution has to be averaged over the stochastic process h{t). We define the 
phase accumulated by m+ as a sum of a regular, 0o and stochastic, (p, parts: 

/ dt'B{t') = Bot+ [ dt'hit') = Mt) + fit) , (4) 
Jo Jo 

and obtain for average value of m_|_ 

(m+(t)) = (e*'^)m+(0) = e^'^°(e^'^)m+(0). 

The Gaussian approximation is based on the assumption that the probability 
distribution of the phase ip is Gaussian: 

1 _ 

PM = ^^=^e (5) 



4>{t) 



As usual, this assumption is to be justified on the basis of the central limit theorem. 
The stochastic phase ip is the integral of the random process 

v{t) ^ b{t) . (6) 

The Bloch vector precesses around the 2;-axis with the angular velocity that has 
random modulation v{t) In the Gaussian approximation the only relevant statistical 
characteristics of v(t) is the correlation function {v(ti)v(t2)) = W{\ti — ^2!) (we 
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assume that v{t) is a stationary random process). This function W{t) vanishes at 
r ^ oo and the scale of this decay is the correlation time, Tc- If the time of 
integration, t, is much longer than Tc the random phase (p{t) is a sum of many 
uncorrelated contributions. The central limit theorem implies that such a sum has 
Gaussian distribution, independently of the details of the process. This is our initial 
understanding: the Gaussian approximation becomes valid as soon as t exceeds the 
correlation time of the noise. Below we will further discuss this conclusion. It follows 
from Eq. ([5]) that 

(e*'^) = j dipp{ip)e'^ = e-^<'^'> (7) 

with 

{if^) = [ dti [ dt2 {v{h)v{t2)) = [ dti [ dt2W{\ti-t2\). (8) 
Jo Jo Jo Jo 

The correlator W{t) = {v(t)v{t + r)) is usually represented by its Fourier transform - 

the power spectrum of the noise, S{ll!): 

S{lj) = - / dtW{t) cos cut. (9) 







TT 

Using Eqs. (jUinD we obtain 



/oo sin^ — 
dcu ^Sicu). (10) 

For large t the identity 

sin^ ax ^ , . 
lim = 0{x) 

implies that (yj^) 27rt5'(0) and thus 

(e*^) = e-*/^^ T2^ = ixS{<d). (11) 

Therefore, the Gaussian approximation leads to the exponential decay of the signal at 
large times, the decrement being given by the noise power at zero frequency. 



Gaussian approximation: Echo experiments 

In the following we will discuss the Bloch vector dynamics in the rotating frame of 
reference where 0o = and m+(t) = e'^'^^^\ The time dependence of = (e**^) 

characterizes decay of the so-called free induction signal [S^ . In order to extract it in 
qubit experiments, one has to average over many repetitions of the same qubit operation. 
Even in setups that allow single shot measurements [35] each repetition gives one of the 
two qubit states as the outcome. Only by averaging over many repeated runs can one see 
the decay of the average as described by the free induction signal. The problem with this 
is that the environment has time to change its state between the repetitions, and thus 
we average not only over the stochastic dynamics of the environment during the time 
evolution of the qubit, but over the initial states of the environment as well. As a result, 
the free induction signal decays even if the environment is too slow to rearrange during 
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the operation time. This is an analog of the inhomogeneous broadening of spectral lines 
in magnetic resonance experiments. This analogy also suggests ways to eliminate the 
suppression of the signal by the dispersion of the initial conditions. One can use well 
known echo technique(see, e. g., [36]) when the system is subject to a short manipulation 
pulse (so-called vr-pulse) with duration ri at the time tu. The duration ti of the pulse 
is chosen to be such that it switches the two states of the qubit. This is equivalent 
to reversing the direction of the Bloch vector and thus effectively reversing the time 
evolution after the pulse as compared with the initial one. As a result, the effect of any 
static field is canceled and decay of the echo signal is determined only by the dynamics 
of the environment. The decay of 2-pulse echo can be expressed as (m^''(2ri2)) |36j . 
where 



\J0 Jti2/ 

Finite correlation time of v{t) again leads to the Gaussian distribution of ilj{t) at large 
enough t with 



This variance can be much smaller than (v?^), Eq. ffTOl) . if S{u!) is singular at ^ 0. 
3. Non-Gaussian decoherence: Spin-fiuctuator model 

We have discussed how to calculate the decoherence of a qubit in the Gaussian 
approximation. The only statistical characteristic of the noise is the time correlation 
function, W{t), or equivalently the power spectrum, S{uj). Now let us consider situations 
when this approximation is not valid. In such situations knowledge of only the noise 
power spectrum is not sufficient: noise sources with identical power spectra can have 
different decohering effect on the qubit. Thus, it is necessary to specify the model for 
the noise source in more detail. In the following we will use a random telegraph process 
as a model noise source. We will first analyze individual random telegraph processes 
and show how 1/f noise appears as a result of averaging over a suitable ensemble of 
these processes. We will describe the decoherence due to one telegraph process, and 
finally extend this discussion to decoherence by averaging over ensembles of telegraph 
processes. 

3.1. Random telegraph processes 

Consider a stochastic function x(^); which at any time takes the values x(^) = il [Zj- 
It is thus suitable for describing a system that can find itself in one of the two stable 
states, 1 and 2, and once in a while makes a switch between them. The switchings are 
assumed to be uncorrelated random events with rates 712 and 721, which in principle can 
be different. Here we will limit ourselves to symmetric telegraph process: 712 = 721 = 7- 




(12) 




(13) 
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The extension to the general case is straightforward, see, e.g., [3Z]. The number k of 
switches that the fluctuator experiences within a time t follows a Poisson distribution 

Pk = ^e-^* . (14) 

The number of switches, k, determines the number of times the function x(t) changes 
its sign contributing (— l)'^ to the correlation function, C{t) = (x('^)x(O))- Therefore 

C(t) = e"^*^(-l)'=^^ =e-2^*, t>0. (15) 

k=0 

The random telegraph process results in fluctuating field v(t) = ^rx(t) applied to the 
qubit. The magnitude of this field, v = |v|, together with the switching rate, 7, 
characterizes the fluctuator. Using (Q we find power spectrum of the noise generated 
by i-th fluctuator: 

Generating 1/f noise by sets of telegraph processes. - Let us consider a combination 
of a few statistically independent telegraph processes. Since the effective fields are 
all parallel to the z-axis we can characterize z-th fluctuator by a coupling strength Vi 
and a switching rate 7, so that the total magnitude of the fluctuating field is given by 
'^iViXiit). The total noise spectrum S{u!) equals to a sum ^j>S'i(a;). If the number of 
fluctuators interacting with the qubit is large the fluctuating field can be written as a 
convolution of Si^uj), Eq. f|T6|) . with the distribution P{v, 7) of the parameters v and 7, 

S(.) = l.jv'ivjd-<Viv,-<)-^^. (17) 

To obtain the l/uj low-frequency behavior of the power spectrum f[T7|) one has to assume 
that the distribution function P{v,'y) behaves as I/7 at small 7, i.e., 

V{v,^)\,^o = V{v)/^. (18) 

In this case, 

POO 

S{uj) = {v^)/uj, {v^)= v^Viv)dv. (19) 

Jo 

It turns out that the distribution (ITSi) . which corresponds to a uniform distribution of 
log 7, follows naturally from a simple model. In this commonly used model the role 
of a bistable fluctuator is played by a particle confined by a double well potential |38j ; 
each well has one state to host the particle, and the tunneling between these states is 
possible. This two level system (TLS) describes a broad class of fluctuators: the particle 
can be a generalized particle moving in a generalized space. Hamiltonian of a TLS with 
tunneling can be written as 

Hi. = i (Ar, + At,) (20) 
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where are the Pauh matrices (in contrast with the Pauh matrices cxj that act in the 
qubit Hilbert space). Each TLS is characterized by two parameters - diagonal sphtting, 
A, and tunnehng matrix element, A = hujQe'^. Here uq is a typical frequency of the 
classical motion of the "particle" inside each of the two wells. The logarithm of the 
switching rate 7 is proportional to the dimensionless tunneling integral A. A natural 
assumption that A is uniformly distributed in the TLS ensemble leads to the distribution 
f fTSjl for the switching rate 7. 

The environment is usually modeled as a thermal bath, which can represent a 
phonon field as well as, e.g., electron-hole pairs in the conducting part of the system. 
Fluctuations in the environment effect the fluctuator through either A or A, Eq. ( l20l) . 
Assuming that the modulations of the diagonal splitting A are most important we can 
describe the interaction of the environment with the qubit as 

T-Cp-env = g'cr^, (21) 

where c is an operator in the Hilbert space of the environment depending on the concrete 
interaction mechanism. It is convenient to diagonalize Tip (l20l) by rotating of the 
fluctuator Hilbert space. Then 

HF = |r,, E = VA2 + A2. 

Keeping the notation Tj for the Pauli matrices representing the fluctuator in the rotated 
basis we write the interaction Hamiltonian fl2Tl) as 



/MAX 

T-tp-env = gc\—T^- —r^ I . (22) 

The factor (A/i?)^ appears in the inter-level transition rate [39| HQ]: 

^{E,A) = {A/E)'^o{E). (23) 

Here the quantity '~fo{E) has the meaning of the maxzma/ relaxation rate for fiuctuators 
with a given energy splitting, E. 

The parameters A and A are usually supposed to be uniformly distributed over 
intervals much larger than those important for low-temperature kinetics [2H1 SI] , so we 
can write their probability distribution as 

V{X,A)=Po (24) 

(leaving aside the problem of normalization, which will be fixed later). It is convenient to 
characterize each fluctuator by two parameters - the energy spacing between the levels, 
E = VA^ + A^, and 9 - determined through the relations A = E cos 6, A = E sin 9 
{9 < 7r/2). The mutual distribution function of these parameters can we written as [12] 



Relaxation rates 7 for the fiuctuators with a given spacing E are distributed according 
to 



or 1 -1/2 

T^{E,1) = — 



27 



7min(^) < 7 < loiE) . (26) 



Decoherence in qubits due to low-frequency noise 



9 



To normalize this distribution one has to cut it off at small relaxation rates at a minimal 
value 7jnin(-E)- The distribution (126|1 has Eq. f|T8l) as its limit at 7 <^ 70 and thus leads to 
1// noise in the frequency domain 7min <^ co" <^ 7o- The coupling, v, between the qubit 
and the fluctuator depends on the parameters E and 9 of the fluctuator, see Sec. 13.31 
for more details. Since only the fiuctuators with E <T are important and both v and 
7o are smooth functions of E one can use the values of v and 70 at -E = T. Below we 
restrict ourselves by this approximation which allows us to perform integration over the 
energy spacings of the fluctuator states. In particular, integrating Eq. (l26l) over E and 
7 we connect Pq with the total number, N^, of thermally excited fiuctuators: 

P -^^ r - In (97) 



The integral expression f[T7j) for the noise power spectrum is valid provided that 
PqT ^ 1, i. e., Mt ^ i^- Substituting Eq. fl271) into Eq. (fT7|) we obtain the estimate: 

S{uj) ~ {v )PqT { . „ (28) 
\ 7o/c^' , > 70 . 

Therefore, the noise in the SF model indeed has the 1// noise power spectrum at low 
enough frequencies. The crossover from uj~^ to uj^'^ behavior at a; ~ 70 manifests the 
cutoff in the distribution fl2B]) at high switching rates. 

3.2. Decoherence by a single random telegraph process 

Master equations. - Having described the properties of random telegraph noise we 
are now ready to discuss how it effects a qubit. Let us turn to Eq. Q with a single 
telegraph process as the noise source. 

We assume that the fluctuator does not feel any feedback from the qubit and thus 
the random telegraph function equals to +1 or -1 with the probability 1/2 regardless 
to the direction of the Bloch vector M. Consider the dynamics of the Bloch vector in 
the rotating frame of reference during a time interval t,t + T, which is small as compared 
to the inverse switching rate I/7, so that the fluctuator changes its state with a small 
probability, 7r ^ 1 while the probability to switch more than once is negligible. 

Let us split the probability to find the angle if at time t, p{if,t), into partial 
probabilities to arrive at the angle (f due to rotation around the vector b and — b, 
respectively: 

p{ip,t)=p+iip,t)+p+iip,t). (29) 

Denoting as V{t) the probability for a fluctuator not to switch during the time interval 
(t, t + r) and the probability to switch once as Q(r) one can express p+{ip,t + T) through 
p+{ip,t) and p-i(p,t): 

p+{^,t + T) =V{r)p+{ip + VT,t) + dti\^Q{ti-t)p4ip-v{t + T-ti),t]y (30) 
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A similar equation can be written for p^{(f, t + t): 

p.{^,t + T)=V{T)p-{^-VT,t)+ dh^Q{ti-t)p+[p + v{t + T-h),t]Y (31) 

When r ^ 7"^ we can indeed neglect multiple switchings, and expand Eq. f|T^ as 
V{t) = = I - -IT, Q{t) = = ^T. Differentiating Eqs. ([30]) and §3) over r we 
arrive at the set of equations: 

p+= - 1P+ + IP- + vd^p+ , 

P- = - + 1P+ - vd^P- , (32) 

d^p± = dp±/d(p. Equations fl32|) can be combined to one second-order differential 
equation for p{(p,t), Eq. ( l29i) . 

p + 2rp = v''dlp, (33) 

which is known as the telegraph equation. We can always choose a;-direction in such 
a way that (p = at t = 0. Thus ^(v^, 0) = S{(p). The second initial condition to 
Eq. fl33|l . p{(p,0) = ±2f(9(pp((y9, 0), follows directly from the integral equation (130|) . The 
sign depends on the initial state of the fluctuator. After averaging over the fiuctuator's 
initial state, p{ip, 0) = 0. 

To evaluate the decoherence for this model we multiply Eq. (l33l) by = e*"^ and 
integrate over ip to show that the mean value of m+ , 



satisfies the equation 

(m+) + 2'j{m+) = -v'^{m+) . (34) 

The initial condition for this equation is m^{0) = 1 since f{0) = 0. The second initial 
condition, m+(0) = 0, follows from the boundary condition p(0, 0) = 0. The solution of 
Eq. (IMl) with these initial conditions, 

{m^) = ^ + 1)^'"* + - 1)^"'"*] ' = ^1-^ , (35) 

describes decoherence of a qubit due to a single random telegraph process given by the 
coupling strength v and the switching rate 7. 

According to (l35l) . the free induction signal demonstrates qualitatively different 
behavior for large and small ratios f/7. At ^ 7 one can consider the qubit as 
a quantum system experiencing beatings between the states with different splittings. 
Bo ± b, the width of these states being 7. In the opposite limiting case, f <^ 7, the 
inter-level splitting is self-averaged to a certain value, the width being f^/27. This 
situation was extensively discussed in connection with the magnetic resonance and is 
known as the motional narrowing of spectral lines [I2] . The two types of behavior will 
be discussed in more detail in Sec. HI 
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Comparison to the Gaussian approximation. - Let us now compare Eq. fl35p with 
the result (fTTj) of the Gaussian approximation. Substituting f|T6|l for the noise power 
spectrum one obtains 



(G) 



2^ 



(36) 



When discussing the Gaussian approximation we argued that it should be valid for 
times longer than the correlation time of the noise, which for the random telegraph 
process is (27)"^. Expanding our Eq. fl35|) at long times we find that it indeed decays 
exponentially (or the oscillations decay exponentially in the case > 7). However the 
rate of the decay is parametrically different from Eq. (!36l) : 



1 



7 — 7 Re 



(37) 



It is easy to check that in the limit w <^ 7 Eq. ( l37l) coincides with the Gaussian result. 
Shown in Fig. [2] are equations fl5Bl) and fl57|) as functions of v for fixed 7 = 1. We see 




Figure 2. Comparison of the decolierence rate for a single random telegrapli 
process and tlie corresponding Gaussian approximation. 

that the Gaussian approximation is only valid in the limit u -C 7. Apparently this 
conclusion is in a contradiction with the previous discussion based on the central limit 
theorem. It looked convincing that according to this theorem p{ip, t) always tends to a 
Gaussian distribution with time- dependent variance provided that the time exceeds the 
correlation time of the noise. 

To resolve this apparent contradiction let us analyze the shape of the distribution 
function, p{ip,t), which follows from Eq. ( 1331) . The solution of this equation with the 
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-e-^* + vt) + - vt)] + ^e"^* [e((^ + vt) - e((/? - vt)] 



2v 



X 



+ Jo (7tVl - 



where /^(x) is the modified Bessel function and 9{x) is the Heaviside step function: 



(3^ 



e{x) 



1, x > , 
0, x<0. 



The distribution function (1551) (shown in Fig. 15] for various t) consists of two delta- 
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Figure 3. The distribution function ([38)) for t 
positive (p is shown, the function is symmetric), v 
delta-functions (not to scale). 



= 1, 5 and 10 (only the part for 
7=1. The arrows represent the 



functions and a central peak. 

The delta-functions represent the finite probability for the fiuctuator to be in 
the same state during time t. As time increases, the weight of the delta-functions 
decreases and the central peak broadens. At large times this peak acquires a Gaussian 
shape. Indeed, at 7t ^ 1 one can use the asymptotic behavior of the Bessel function 
Iv{z) (27r)~^/^e^, as z oo. For (p <^ vt we can also expand — {(p/vt)"^ and 
convince ourselves that the central peak in accordance with the central limit theorem 
is indeed described by the Gaussian distribution ([5]) with the variance {ip"^) = v'^t/'-f. If 
the qubit-fiuctuator coupling is weak, ?; <^ 7, this Gaussian part of p{ip, t) dominates 
the average (e*''') and the Gaussian approximation is valid. Contrarily, for the strong 
coupling, f > 7, the average is dominated by the delta functions at the ends of the 
distribution, and the decoherence demonstrates a pronounced non-Gaussian behavior, 
even at t > 7^^. 

I Note that in Ref. [30] the term with Iq was missing. 
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Unfortunately it is impossible to measure the distribution of if experimentally. 
Indeed, each experimental shot corresponds to a particular realization of the random 
external noise and does not yield a particular value of (p. The reason is in the difference 
between the qubit that can can be viewed as a pseudospin 1/2 and a classical Bloch 
vector M. According to Eq. ([2]) the components Mx, My, Mz of M are connected with 
the mean component of the final state of the pseudospin. Therefore, to measure the 
value of the phase ip (argument of m_|_) that corresponds to a given realization of the 
noise one would have to repeat the experimental shot with the same realization of the 
noise many times. This is impossible because each time the realization of noise is 
different. Therefore, the only observable in the decoherence experiments is the average 
(e*'^). There is no way to extract more information about the distribution p(v5, t) from 
any experiment with a single qubit. 



3.3. Averaging over ensembles of random telegraph processes 

General expressions. - As we have seen, a set of fluctuators characterized by the 
distribution function ( 1261) of relaxation rates provides a realistic model for 1/f noise. 
Thus, in order to study the decoherence by such noise it is natural to sum contributions 
of many fluctuators. To perform this procedure we assume that dynamics of different 
fluctuators are not correlated, i.e., (Xj(^)Xi(^')) = 5jje~^'^'l*~*'l. Under this assumption 
the average value of the complex momentum, (m+) is just a product of the partial 
averages, 

(m+(t)) = J](m+,(t))=e^»''^<"^+»«>. 

i 

Since the logarithm of a product is a self-averaging quantity, it is natural to approximate 
the sum of logarithms, ln(m+(t))j, by its average value, —]Cm{t), where 



/C„(t) = -5^1n(m+,(t)). (39) 

i 

Here bar denotes the average over both the coupling constants, v, of the fluctuators 
and their transition rates, 7. If the number At of thermally excited fluctuators is large 
we can replace the sum ^jln(m+i) by A/'Tln(m+). Furthermore, one can employ the 
Holtsmark procedure [13], i. e., to replace ln(m_|_) by (m^) — 1, assuming that each of 
(m_|_i) is close to 1. Thus, ICmif) is approximately equal 



lCm{t) ^ Mt{1 - {m+)) = J dvd-fV{v,-f)[l- {m+{v,^\t))] . (40) 

Here (m_|_(w, 7|t)) depends on the parameters v and 7 according to Eq. (!35|l . The average 
free induction signal is then (m_|_(t)) = e~''^'"^*\ 

Analysis of the echo signal is rather similar: one has to replace (m_|_(t)) taken from 
Eq. (ESD by Hg 



-27ri2 



2?;2' 

(/i + l)e2^^"^^ -ifi- l)e-2^^^i2 - -- 



(41) 



Decoherence in qubits due to low-frequency noise 



14 



To evaluate the time dependence of either the free induction or the echo signal 
one has to specify the distribution of the coupling constants v. Let us consider each 
fluctuator as a two-level tunneling system characterized by the Hamiltonian 

^? = ^(A.r«+A.r«) (42) 

where t*-*-* is the set of Pauli matrices describing i-th. two-state fluctuator. The energy 
splitting for each fluctuator is Ei = Aj + Kf. The variation of the energy spacing 
between the states of the qubit can be cast into the effective Hamiltonian, which (after 
a rotation similar to that in Eq. fl2^ ) acquires the form 

HqF = ^Vi f^zT^'^ , Vi = g{ri)A{ni) cos 9i . (43) 

i 

Here 6i = arctan(Aj/Aj), rij is the direction of elastic or electric dipole moment of i-th 
fluctuator, and is the distance between the qubit and i-th fluctuator. Note that in 
Eq. (H3l) we neglected the term proportional to azTx- This can be justified as long as 
the fluctuator is considered to be a classical system. The functions A{ni) and g{ri) are 
not universal. 

The coupling constants, fj, determined by Eq. (143|) . contain cos6'j and thus are 
statistically correlated with 6i. It is convenient to introduce an uncorrelated random 
coupling parameter, Ui clS 

Ui = g{ri)A{ni) , Vi = Ui cos di . (44) 

It is safe to assume that direction, rij, of a fluctuator is correlated neither with its 
distance from the qubit, rj, nor with the tunneling parameter represented by the variable 
6i. This assumption allows us to replace ^(n) by its angle average, A = (|yl(n)|)n. 




Figure 4. Sketch of localized charges near an electrode. Induced miage charges create 
local dipoles which interact with the qubit. 

Simple model. - The coupling parameter, g, decays with distance, r, between the qubit 
and the fluctuator. Usually this decay is algebraic: g{r) oc g/r'^. We will distinguish 
between two cases: (i) the fluctuators are distributed in three dimensional space (d=3) 
and (ii) the fluctuators are located in the vicinity of a two dimensional manifold, e.g.. 
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in the vicinity of the interface between an insulator and a metal (d=2). Using the 
distribution (126|) of the relaxation rates one can express V{u, 9) as 

nM) = (.co.9rjI^^. (45) 

Here is a rf-dependent dimensionless constant while is a typical distance between 
the fluctuators with Ei < T. In the following we will for simplicity assume that 

^min < < r 

max ) 

(46) 

where rmin (^'max) are distances between the qubit and the closest (most remote) 
fluctuator. Under this condition rj cx T^^'^ is the typical constant of the qubit-fluctuator 
coupling. As soon as the inequality (l46l) is violated the decoherence starts to depend 
explicitly on either rmin or rmax, i- e., becomes sensitive to particular mesoscopic details 
of the device. 

Let us first consider the case when ci = 6, as it is for charged traps located near 
the gate electrode (were b = d = 2 due to the dipole nature of the field produced by a 
charge and its induced image, see Fig. |4]). In this case one can cast Eq. (140|) into the 
form 

]Cf{t) =7] — d9 tan ^ [1 - f{u cos 9, 70 sin^ 9\t)\ . (47) 

Here /(f,7|t) is equal either to (m+(t>,7|t)) or to (m^''(t>, 7|t)), depending on the 
manipulation protocol. Equation (1471) together with Eqs. (135|) and (14T1) allows one 
to analyze different limiting cases. 

In the case of two-pulse echo it is instructive to look at the asymptotic behaviors 
of }Cf{t). From Eq. (gl]) it follows that 

{t^7osin6'(Mcos6')^ , t < (70 sin 61)-^ 
t^{ucos9f , (7osin^)-i < t < (^43) 

tucos9 , u^^ <^ t. 

Splitting the regions of integration over u and 9 according the domains (HHj) of different 
asymptotic behavior one obtains 

/Cm(2ri2) ~ 77x12 min{7ori2, 1} . (49) 

The dephasing time (defined for non-exponential decay as the time when /C ~ 1) 
for the two-pulse echo decay is thus equal to 

r^ = max{r/-\(r/7o)-i/2} . (50) 

The result for 70T12 <^ 1 has a clear physical meaning [12]: the decoherence occurs 
only provided that at least one of the fluctuators flips. Each flip provides a contribution 
~ Tjt to the phase, while 70T12 ^ 1 is a probability for a flip during the observation 
time. The result for 7ori2 ^ 1 is less intuitive since in this region the dephasing is 
non-Markovian, see Il2] for more detail. 
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It is important that at large observation times, T12 ^ % , the decoherence is 
dominated by few optimal fluctuators . The distance ropt(T), between the optimal 
fluctuators and the qubit is determined by the condition 

^^(ropt) ~ 7o(T) . (51) 

Derivation of Eq. (ISTl requires a rather tedious analysis of the expansion (148|) and the 
integration over u and 6, see |^ [25] for more detail. This estimate emerges naturally 
from the behavior of the decoherence in the limiting cases f ^ 7 and v <^ 7. For strong 
coupling the fluctuators are slow and the qubit's behavior is determined by quantum 
beatings between the states with E±v. Accordingly, the decoherence rate is of he order 
of 7. In the opposite case, as we already discussed, the decoherence rate is oc t'^/7. 
Matching these two limiting cases one arrives at the estimate (l5Tll . 

What ii d ^ hi If the coupling decays as l/r'' and the fluctuators are distributed 
in a d-dimensional space, then r'^~^ dr — > V{v) oc As a result, at d < b the 

decoherence is controlled by optimal fluctuators located at the distance ropt provided 
that they exist. At d > b the decoherence at large time is determined by most remote 
fluctuators with v = iVn- li d < b, but the closest fluctuator has t'max ^ 70) then 
it is the quantity fmax that determines the decoherence. In the last two cases /C(t) is 
proportional to t^, and one can apply the results of [22], substituting for v either v^in 

or f^ax' 

Since ropt depends on the temperature one can expect crossovers between different 
regimes as a function of temperature. A similar mesoscopic behavior of the decoherence 
rate was discussed for a microwave-irradiated Andreev interferometer |44j . 

Note that the result ( l49i) for the long-range interaction cannot be reproduced by the 
Gaussian approximation since in the latter case the decoherence would be determined 
by the nearest neighbors of the qubit. That can be seen from explicit expression for 
the phase accumulation following from the Gaussian model, see Eq. (l52l) below. Being 
substituted in Eq. ( 1471) this expression leads to divergence of the integral over u at its 
upper limit that physically means dominating role of nearest neighbors of the qubit. At 
the same time, the spin-fluctuator model implies that the most important contributions 
are given by the fluctuators with v{r) ~ 70. 

The above procedure still leaves unanswered a delicate question: Can contributions 
of several fluctuators be described by averages over the fluctuators' parameters? The 
situation with a qubit interacting with environment in fact differs from that of a resonant 
two-level system in spin or phonon echo experiments. In the first case the experiment is 
conducted using a single qubit surrounded by a set of fluctuators with fixed locations, 
while in the second case many resonant TLSs participate the absorption. Consequently, 
one can assume that each TLS has its own environment and calculate the properties 
averaged over positions and transition rates of the surrounding fluctuators. How many 
surrounding fluctuators one needs to replace the set of fluctuators with fixed locations 
(and transition rates) by an averaged fluctuating medium? This issue was studied 
numerically in [25]. The analysis showed that one needs really many (> 100) fluctuators 
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4. Relevance to experiments 

In this section we will briefly review some experiments were we believe that the 
theory of non-Gaussian noise is relevant. While there are some indications of non- 
Gaussian behavior, to our knowledge, the existing experiments are not conclusive 
enough. However, we should emphasize that observation of non-Gaussian effects was 
not among the goals of thee experiments. The main goal was to achieve the longest 
possible decoherence times. Probably as progress is made the situations where only 
one or a few independent noise sources are important would become rather usual than 
exceptional making the non-Gaussian effects more pronounced. We also believe that 
experimental studies of the non-Gaussian effects would provide useful information on 
the environmental degrees of freedom. Therefore devices showing pronounced non- 
Gaussian behavior, while not necessarily being the optimal qubits, may serve as useful 
research tools. 

4.I. Plateaus in the echo signal 

Using Eqs. ( fTTI) and ( l36l) we have calculated the decay of the free induction signal due 
to a single fluctuator in the Gaussian approximation at large times. Let us return to 
Eq. (fT3!) and evaluate the same for the echo signal, but at arbitrary times. The result is 



with (m^ ) = e"^*^^^/^ as before. This should be compared with the exact result, 
Eq. (1411) . In Fig. \5\ both results are shown for both a weakly (f/7 = 0.8) and a 
strongly {v/'-f = 10) coupled fluctuator. As has been discussed before, the Gaussian 




^ [47ri2 - 3 + 4e"2^"^^ 




(52) 




0, 







0.5 



1.5 



2 



ft 



Figure 5. Echo signal for different values of the ratio u/7 (shown near the curves), 
Eq. (|¥T|l . Dashed lines - calculations along the Gaussian approximation, Eq. ([5^ . 
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approximation is accurate for f < 7, while at f > 7 the two results are qualitatively 
different. In particular, the plateaus in the time dependence of the echo signal shown 
in Fig. [5] are beyond the Gaussian approximation. We believe that such a plateau was 
experimentally observed in Ref. [15] (see Fig. 3 there). In the limit w S> 7, \/'y/T, 
equation (HTl) acquires a simple form: 

(mfh = fl + - sin2t;rl . (53) 

L V 

According to Eq. fl53l) . the plateau-like features {d{m!'^) / dr ^ 0) occur at vt ~ kn 
and their heights (m^^) g-27rfc7/f exponentially decay with the number k. Measuring 
experimentally the position and the height of the first plateau, one can determine both 
the fluctuator coupling strength v and its switching rate 7. For example, the echo 
signal measured in Ref. |15] shows a plateau-like feature at T12 = 3.5 ns at the height 
(m^^) = 0.3, which yields v ^ 143 MHz, and 7 f« 27 MHz. If the fluctuator is a 
charge trap near the gates producing a dipole electric field its coupling strength is 
V = e^(a • r)/r^. Using the gate-qubit distance r ^ 0.5 /im, we obtain a reasonable 
estimate for the tunneling distance between the charge trap and the gate, a ~ 20 A. A 
more extensive discussion of this is found in Ref. [281 ■ 



4.2. Flux qubit 

In Sec. 13.31 we discussed a model with a broad distribution of coupling strengths, v. 
This is appropriate in a situation were the noise sources are distributed uniformly in 
space and act on the qubit via a long range (power law) force. In this section we want 
to apply the spin-fiuctuator model to experiments on flux qubits [^ HTJ. Since the 
microscopic source of flux noise is not clarified (see, e.g., [H] and references therein), 
it is not clear what would be the most reasonable distribution of v. We have adopted 
the simplest model, where the coupling parameters v are narrowly distributed around 
some characteristic value v. In other words, v and 7 are supposed to be uncorrelated, 
P(w,7) =P,(^;)P^(7), and 

P T 

^7(7) = -^<^(7o-7)<^(7-7min), V^{v) = 5{v -v). (54) 

7 

Using Eq. (JSlj) and the expressions for (m^^) for either the Gaussian [Eq. (l52l) ] or the 
spin-fiuctuator [Eq. fHT]) ] models we obtain the quantity /C defined in Eq. fl39l) . In the 
Gaussian model: 

ir fn \ nr>rr-2 2 1 270X12/3 , 7ori2 < 1 
/Cg 2ri2 = 2PoTv r^a x < , ^ ^1 

In 2 , 7ori2 > 1 . 

We see that at long times, 7ori2 ^ 1, we have a quadratic dependence on time, which 
manifests a Gaussian decay of the echo signal. At short times, 70X12 <^ 1, the result is 
multiplied by an additional factor 7or, which is the probability for a single flip of the 
fastest fiuctuators. 

In the spin-fiuctuator model there are two limiting cases: 
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(i) When v 



^sf(2ri2) 



PoT(21n2)t;V2 
PoTavTi2 



Tl2 < 7(7\ 
V"^ < Tl2- 



(56) 



where a ~ 3. At small times ri2 <^ we arrive at the same result as in the 
Gaussian approach, Eq. ( l55l) . However, at large times, ri2 ^ -w"^, the exact 
result dramatically differs from the prediction of the Gaussian approximation. To 
understand the origin of the non-Gaussian behavior notice that for ^7(7) oc I/7, 
the decoherence is dominated by the fluctuators with v. Indeed, very "slow" 
fluctuators produce slow varying fields, which are effectively refocused in course of 
the echo experiment. As to the "too fast" fluctuators, their influence is reduced 
due to the effect of motional narrowing. We have already learned that only the 
fluctuators with <^ 7 produce Gaussian noise. Consequently, the noise in this 
case is essentially non-Gaussian. Only at times T12 ^ v"^-, which are too short for 
these most important fluctuators to switch, the decoherence is dominated by the 
faster fluctuators contribute, and the Gaussian approximation turns to be valid. 

(ii) When w ^ 70 we find 



In this case all fluctuators are strongly coupled to the qubit. Therefore, the long- 
time decoherence is essentially non-Gaussian. 

We thus conclude that it is the long-time decoherence that is most sensitive to the 
particular model of the noise. Unfortunately, at long times the signal usually is weak 
and obscured by noise. One of the possible ways to experimentally identify the non- 
Gaussian behavior is based on the fact that the typical fluctuator strength, -u, enters the 
the expressions for the Gaussian and non-Gaussian decay in different ways. In particular, 
according to the Gaussian approximation, }Cm{t) oc v'^ at all times. Contrarily, according 
to the SF model [Eqs. ( l56l) . (1571) ] the powers of v are different for different times. By 
fitting to the spin-fluctuator model this parameter can be extracted, and from this we 
can infer the flux change corresponding to the flip of a single fluctuator. We have 
performed such fits to the experimental data of [16], the details are given in [17], and 
one example is reproduced in Fig. [61 As can be seen, it is difficult to determine on the 
basis of these experiments whether non-Gaussian effects are important, or the Gaussian 
approximation is sufficient, both models fit the data equally well. However, if we choose 
to fit the spin-fluctuator model and fit not only to one experimental curve, but to the 
whole set of curves for different working points of the qubit, we find that Eq. fl56l) for the 
case y <^ 7o provides the better overall fit. The change of flux due to a single fluctuator 
flip is < 10~^$0) where $0 is the flux quantum. 




T12 < V 
T12 > V 



-1 



(57) 
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time (|is) 

Figure 6. Fit of the experimental data from [46] to the Gaussian model, and to the 
spin-fluctuator model with v <^ [Eq. ((56)) ] or 5 ^ 70 [Eq. ([57)) ] . 

5. Microscopic sources of telegraph noise 

In this section we briefly discuss of microscopic noise sources that can produce classical 
telegraph noise. 

Charge noise - The obvious source of RTP charge noise is a charge, which jumps 
between two different locations in space. Less clear is where these charges are actually 
located and what are the two states. The flrst attempt of constructing such a model 
in relation to qubit decoherence appeared in [22], where electrons tunneling between a 
localized state in the insulator and a metallic gate was studied. This model has been 
further studied in [19], [50] . Later, experimental results [1] indicated a linear dependence 
of the relaxation rate on the energy splitting of the two qubit states. One also has to take 
into account that in the experimental setup there is no normal metal in the vicinity of the 
qubit: all gates and leads should be in the superconducting state at the temperatures of 
experiment. These two facts suggest that the model [22] is irrelevant for the decoherence 
in charge qubits [15] and favored a model with superconducting electrodes [51]. In this 
model, the two electrons of a Cooper pair are split and tunnel separately to some 
localized states in the insulator (see Fig. [7] for an illustration of this (Model III) and 
other models). A constant density of these localized states gives a linearly increasing 
density of occupied pairs, in agreement with experiments [1]. This model was criticized 
[52] because it required an unreasonably high concentration of localized states, and a 
more elaborate model was proposed. However, it was shown [IHl [50] that allowance for 
quantum effects of hybridization between the electronic states localized at the traps and 
extended states in the electrodes relaxes the above requirement. At present it seems 
that no solid conclusions can be drawn based on the available experiments. 
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Model I Model II Model III 



Figure 7. Three possible models for the fluctuating charges: Model I, electrons 
jumping between a localized state and a normal metal, as discussed in [22l |49l [50] . 
Model II, electrons jumping between localized states. Model III, electrons jumping 
between localized states and a superconductor, as discussed in |51j . 

Noise of critical current - The microscopic mechanism and the source of the 
fluctuations on the critical current in a Josephson junction is a long-standing open 
problem. These fluctuations were initially attributed to the charges tunneling or 
hopping between different localized states in the barrier forming glass-like TLSs. 
However, a more detailed comparison with experiment revealed an important problem 

- the experimentally observed noise spectrum [53] was proportional to that is 
incompatible with the assumption of constant TLS density of states. The interest in 
the critical current fluctuations was recently renewed because of their importance for 
superconducting qubits. The new experiments [51] on fluctuations in small Al junctions 

- similar to those used in several types of qubits [55] - in normal state brought a new 
puzzle. It turned out that the temperature dependence of the noise power spectrum 
in the normal state is linear, and the noise power is much less than that reported 
for large superconducting contacts. A plausible explanation of such behaviors is given 
in [56] where it is suggested that the critical current noise is due to electron trapping in 
shallow subgap states that might be formed at the superconductor-insulator boundary. 
This mechanism is similar to that suggested earlier for the charge noise |52j . 

Flux noise - Studies of the flux noise in superconducting structures have a long history. 
As early as in the 1980s it was demonstrated it is the flux and not the critical current 
noise that limits the sensitivity of dc SQUIDs [57]. The interest in this problem was 
recently renewed when it was realized that flux noise can limit the coherence in flux and 
phase superconducting qubits [^ [58] . Two recent models for fluctuators producing 
low-frequency noise were suggested. The first one [59] attributes the flux noise to the 
electron hopping between traps in which their spins have fixed, random orientations. The 
second one [60] proposes that electrons flip their spins due to interaction with tunneling 
TLSs and phonons. These models were recently criticized in Ref. [61] where it was 
stated that it is hard to justify the assumptions behind both models. In that paper, a 
novel mechanism - based on spin diffusion along the surface of a superconductor - was 
suggested. This model seems to agree with recent experiments [62] on measurements of 
the 1// noise. It remains to be understood whether the surface spin fluctuations lead 
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to pronounced non-Gaussian behavior in the qubit decoherence. 
6. Conclusions 

In this review we have discussed the spin-fluctuator model for qubit dephasing. This 
model provides a simple, solvable, yet in many situations realistic model of the qubit 
environment. In particular, we show how to apply the model in situations where the 
noise shows 1/f behavior. The model shows pronounced non-Gaussian behavior, and 
thus serves as an example of possible deviations from the Gaussian approximation, as 
well as shedding light on the limitations of the Gaussian approximation. 

The main results obtained in studying this model can be summarized as: 

• A single fluctuator, characterized by the coupling strength v and the switching rate 
7 can be classified as weak {v < 7) or strong (f > 7). For weak fluctuators the spin- 
fiuctuator model reproduces the Gaussian result, whereas for strong fluctuators it 
shows non-Gaussian behavior 

• This non-Gaussian behavior persists even in the limit of long time and for an 
arbitrary number of independent fluctuators, as long as each fluctuator is strong. 
This can be understood as a consequence of the 5-functions at the extreme of the 
phase distribution function. 

• In the non-Gaussian case, the time correlation function of the noise is not sufficient 
to determine the qubit dephasing and a more detailed model must be specified. 

• In the non-Gaussian case, individual fluctuators leave signatures in the measured 
signal, e. g., plateaus in the echo signal, that can be used to identify the fluctuator 
parameters. 

We should emphasize that a strong fluctuator, giving non-Gaussian effects, does 
not imply strong decoherence and therefore a bad qubit performance. It indicates only 
that the coupling is strong relative to the switching rate. However, independently of 
the performance of the device as a working qubit, i. e., on whether it achieves a long 
dephasing time, we believe that searching for signatures of non-Gaussian behavior can 
provide valuable information on the nature of the noise source. Given the present 
uncertain state of understanding of the microscopic sources of noise in most solid 
state qubit devices, this seems an important undertaking, and analyzing experiments 
according to the formulas we have presented can be a useful tool in this process. 
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